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ABSTRACT 

We present the Kepler photometric light-variation analysis of the late-type double- 
lined binary system V568 Lyr that is in the field of the high metallicity old open 
cluster NGC 6791. The radial velocity and the high-quality short-cadence light curve 
of the system are analysed simultaneously. The masses, radii and luminosities of the 
component stars are Mi = 1.0886 ± 0.0031 M©, M 2 = 0.8292 ± 0.0026 M©, i?i = 
1.4203 ± 0.0058 RQ, R 2 = 0.7997 ± 0.0015 RQ, L x = 1.85 ± 0.15 LQ, L 2 = 0.292 ± 
0.018 LQ and their separation is a — 31.060 ± 0.002 7?©. The distance to NGC 6791 
is determined to be 4.260 ± 0.290 kpc by analysis of this binary system. We fit the 
components of this wcll-detached binary system with evolution models made with the 
Cambridge STARS and ev(twin) codes to test low-mass binary star evolution. We find 
a good fit with a metallicity of Z = 0.04 and an age of 7.704 Gyr. The standard tidal 
dissipation, included in EV(twin) is insufficient to arrive at the observed circular orbit 
unless it formed rather circular to begin with. 

Key words: binaries: close - binaries: eclipsing - stars:fundamental parameters - 
stars: individual: V568 Lyr -stars: low-mass - open clusters and associations: individ¬ 
ual: NGC 6791 


1 INTRODUCTION 

Accurate photometric and spectroscopic observations of bi¬ 
nary systems provide well-determined parameters of the 
components. These parameters, in turn, allow us to test 
stellar evolution models. Light variations of a large number 
of variable stars have been obtained by successful observa¬ 
tions of the Kepler space telescope. The open stellar cluster 
NGC 6791 is one of the four open clusters that fall in the 
observed field of the Kepler project (Koch et al. 2010). It is 
the oldest open clusters in the Galaxy and contains many 
close and wide binary systems. However, only a few of them 
fall within Kepler's view. Studies of binary systems that are 
members of a cluster allow us to estimate both the age of 
the host cluster and its distance. 

Galactic clusters provide information about the evolu¬ 
tion of the host galaxy. In this regard, old open clusters are 
crucial and NGC 6791, Be 17, M67, NGC 188, etc. (Phelps 
1997; Bedin et al. 2008a; Meibom et al. 2009; Yakut et al. 
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2009; Basu et al. 2011; Brogaard et al. 2011) are important 
old open clusters. Binary systems that are member of clus¬ 
ters with different chemical composition are good targets to 
test how stellar evolution varies with composition, when we 
can determine their physical parameters and construct ap¬ 
propriate evolutionary models. The system V568 Lyr resides 
at the main-sequence turn-off in the colour-magnitude di¬ 
agram (CMD) of the cluster and this turn-off point gives 
us information on the age of the cluster. Hence, the study 
of this particular system gives substantial information on 
both the stars and the cluster in which they reside. This is a 
crucial system, particularly given its detached configuration 
and accurate radial velocity curve and photometry. 

Boesgaard, Jensen, & Deliyannis (2009) analysed the 
chemical composition of the stars at the turn-off of 
NGC 6791 with Keck + HIRES spectra and estimated 
the metallicity to be more than twice that of the Sun. 
Recently, Boesgaard, Lum, & Deliyannis (2014) repeated 
the study and similarly estimated [Fe/H] = 0.30 ± 0.02. 
Previous studies on the metallicity of the cluster give 
the [Fe/H] ratio as 0.30 — 0.40 (Peterson & Green 1998; 
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Origlia et al. 2006). Bragaglia et al. (2014) analysed archival 
Keck/HIRES spectra of main-sequence turn-off and evolved 
stars of NGC 6791. They obtained C, N, O and Na abun¬ 
dances. There are various methods to determine the age 
of the cluster but determination of the chemical composi¬ 
tion is always an important factor. Brogaard et al. (2011, 
2012) used binary system properties of V568 Lyr to esti¬ 
mate the age of the cluster to be 8.3 ± 0.3 Gyr consistent 
with other estimates of 8 — 10 Gyr (Brogaard et al. 2012; 
Boesgaard, Lum, & Deliyannis 2014). There have been a 
number of photometric studies of NGC 6791 and some 
of its stars. Montgomery et al. (1994) observed the clus¬ 
ter with the UBVI filters of the 0.9 m KPNO telescope. 
Stetson, Bruntt, & Grundahl (2003) obtained detailed BVR 
observations with the 2.5 m Nordic Optical Telescope and 
obtained data for colour-magnitude and colour-colour di¬ 
agrams. Platais et al. (2011) present comprehensive cluster 
membership and g'r' photometry of the cluster. Using Ke¬ 
pler photometry, Corsaro et al. (2012) studied solar-like os¬ 
cillations of red giants of the cluster. Photometric studies of 
some of the variable stars in NGC 6791 have also been made 
by Brogaard et al. (2011), Kaluzny & Rucinski (1993) and 
de Marchi et al. (2007). 

V568 Lyr is a relatively long-period system so 
a full and accurate light curve of the system could 
not be obtained. Kaluzny & Rucinski (1993) and 
Rucinski, Kaluzny, & Hilditch (1996) determined its orbital 
period. Later Grundahl et al. (2008) and Brogaard et al. 
(2011) performed light and radial velocity analy¬ 
ses of ground-based data and determined masses of 
1.0868 ± 0.0039M© and 0.8276 ± 0.0022M©, radii of 
1.397 ± O.O13R0 and 0.7813 ± 0.0053R© for the primary 
and secondary components. 

The held of NGC 6791 is rich in binary systems, most 
of which are members of the cluster. The systems V519 Lyr, 
V522 Lyr and V564 Lyr are interacting close binary systems. 
There is no published spectroscopic study of these systems. 
V568 Lyr is, on the other hand, a well-detached system and 
has excellent radial velocity curves (Grundahl et al. 2008; 
Brogaard et al. 2011). This, together with the new Kepler 
light curve, makes it important to determine accurate phys¬ 
ical parameters of its components and so provide better dis¬ 
tance and age estimates. Here, owing to their high sensi¬ 
tivity and quality, we use public Kepler observations sup¬ 
plemented with the radial velocity and abundance measure¬ 
ments of Grundahl et al. (2008) and Brogaard et al. (2011). 
In the next section, we discuss the public Kepler data and 
the light curve and radial velocity data of V568 Lyr are 
solved simultaneously. In the third section, the parameters 
of the system and the distance to the cluster are determined. 
In the fourth section, non-conservative binary evolutionary 
models are constructed with ev(twin) and then these are 
compared with those of earlier studies. 


2 ANALYSIS OF THE KEPLER DATA AND 
LIGHT CURVE MODELLING 

The Kepler space telescope has observed nearly 150000 stars 
with different properties, such as planetary components, pul¬ 
sating stars, etc. The accuracy of the Kepler telescope is very 
high. Part of the NGC 6791 field is in the Kepler field. The 


Table 1. Basic properties of V568 Lyr (KIC 2437452). Magni¬ 
tudes, temperature and E(B — V) are from Montgomery et al. 
(1994) and Brogaard et al. (2011) and other parameters from the 
Kepler Input Catalogue (KIC) and SIMBAD. 


Kepler ID 

2437452 

2MASS ID 

19205427+3745347 

<*2000 

19:20:54.28 

<52000 

37:45:34.7 

2MASS J 

15.513 mag 

2MASS H 

15.057 mag 

2MASS Ks 

15.076 mag 

Kp(Kepler) 

16.981 mag 

V 12 

17.54 mag 

Teff/K 

5645 

E(B-V) 

0.16 

Spectral type 

G5V + K3V 



B-V 


Figure 1. The CMD of the open cluster NGC 6791. Data are 
from Stetson, Bruntt, & Grundahl (2003). The large filled circle 
indicates V568 Lyr, located at the main-sequence turn-off in the 
cluster CMD. 

binary system V568 Lyr is shown on the CMD of NGC 6791 
in Fig. 1. Because it lies right at the main-sequence turn¬ 
off, its analysis is very important for the cluster as a whole. 

V568 Lyr was observed during quarters Q1-Q17 with 

long cadence (exposure time of about 30 min) and one quar¬ 
ter (Q10) in short cadence (lmin). The observed parame¬ 
ters are summarized in Table 1. For the light-curve analysis, 
we selected the short-cadence Q10 public data that cover 
30.11 d continuously. This covers more than twice the orbital 
period. The short cadence Kepler observations of V568 Lyr 
are shown in Fig. 2. The duration of the primary and sec¬ 
ondary minima are almost same and are about 7.8 h. The 
secondary minima show a transit, and its duration is about 
2.8 h. These properties of the light curve allow us to deter¬ 
mine rather accurate relative radii of the components. 

As well as the Kepler light curve, accurate radial ve- 
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Figure 2. Short-cadence Kepler observations of V568 Lyr. 


locity observations of the system exist in the literature. 
Grundahl et al. (2008) and Brogaard et al. (2011) solved 
radial velocities with part of the light curve of V568 Lyr 
obtained with ground-based telescopes. Here, we use the 
highly accurate and full light curve obtained by Kepler. The 
system analysis is made simultaneously with Kepler pho¬ 
tometric observations and the radial velocity curves. We 
expect to obtain more accurate orbital and physical pa¬ 
rameters. The Kepler are solved with the phoebe code 
(Prsa & Zwitter 2005), based on the Wilson & Devinney 
(1971) code. Initial parameters for the light-curve analysis 
are taken from Grundahl et al. (2008) and Brogaard et al. 
(2011). The gravity-darkening coefficients g\ and g 2 are 
obtained from Lucy (1967). Albedos A\ and A 2 are from 
Rucinski (1969) and the limb-darkening coefficients from 
Claret & Bloemen (2011). The orbital inclination i, mass 
ratio q, temperature of the secondary T 2 , surface potentials 
fli and fl 2 , luminosity L\, phase shift <f), the time of mini¬ 
mum light To and the orbital period P are treated as free 
parameters. Previous observations have suggested a third 
body in the system so, during the light-curve analysis, the 
contribution of a third body to the total light parameter 
I 3 is added as a free parameter. With a detailed solution, 
the orbital, geometric and radiative parameters we obtain 
are given in Table 2. The results obtained from observations 
(dots) and from our light-curve model (solid line) are shown 
in Fig. 3a. To show the accuracy of the solution the residuals 
and minima are expanded in Figs 3b-d. In Fig. 3e, the radial 
velocity curves of the components and those calculated with 
simultaneous solutions are shown. 


3 PHYSICAL PARAMETERS OF THE 
SYSTEM 

Light curve analysis of V568 Lyr gives radii, as fractions of 
the separation, Ri/a = 0.04568 and Ri/a = 0.02572. These 
small ratios mean that the components closely preserve their 
spherical symmetry. Because of the high sensitivity of the 
Kepler data, the light curve of the system gives us the op¬ 
portunity to determine much more accurately its inclination 
and especially the radii of the components relative to pre- 


Table 2. Photometric and spectroscopic elements of V568Lyr 
and their formal la errors. See text for details. 


T 0 as JD/d - 2400000 

55741.71981(46) 

P/d 

14.47672(50) 

i/° 

89.05(4) 


21.71(5) 

^2 

32.87(3) 

q = m 2 /mi 

0.7616(1) 

Ti/K 

5645 

Tit K 

4734(80) 

Luminosity ratios: 


^1 Atotal 

74.5% 

^ 2 /^total 

8.5% 

^3 Atotal 

17.0% 

Fractional radius of primary ( Ri/a ) 

0.04568(11) 

Fractional radius of secondary {R 2 / 0 ) 

0.02572(10) 

K i/km s —1 

46.92(9) 

A2/kms -1 

61.60(12) 

System velocity V^/kms -1 

-46.6(1) 

e 

0.000(0.001) 

a/RQ 

31.060(2) 


vious studies. Third light is detected at the level of 17 per 
cent. The reddening is computed, from E(B — V) in Table 
1, to be Ay = 0.49. Radial velocity curves yield the velocity 
amplitudes of the components as listed in Table 2. 

When we leave the eccentricity as a free parameter dur¬ 
ing the solution it tends to zero. Using the orbital parameters 
(Table 2) we obtained the physical parameters of V568 Lyr 
and list them in Table 3. While solving, the temperature 
of the Sun is taken to be 5772 K and its bolometric mag¬ 
nitude to be 4.755 mag (Pecaut & Mamajek (2013)). The 
temperatures of the components with the Flower (1996) ta¬ 
bles give the bolometric corrections as BCi = —0.10 mag, 
and BC 2 = —0.44 mag. Using these BC values for V band, 
we then obtain the absolute magnitudes of the components. 

These new estimates of the system parameters, with 
the magnitudes and reddening, reveal the distance of the 
system to be 4.26 ± 0.29 kpc. The masses we find are consis¬ 
tent with those found by Brogaard et al. (2011). The radius 
of the primary star, on the other hand, is estimated to be 
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0.75 1 _l_l_l_I_l_l_l_I_l_l_l_I_l_l_L 


0.25 0.50 0.75 1.00 1.25 


Table 3. Absolute parameters of V568 Lyr. Standard errors of 
1<7 in the last digits are given in parentheses. 



Primary 

Secondary 

Masses M/M0 

1.0886(31) 

0.8292(26) 

Radii R/ R© 

1.4203(58) 

0.7997(15) 

Effective temperatures T e ff/K 

5 645(95) 

4 734(80) 

Luminosities L/L© 

1.846(15) 

0.292(18) 

Surface gravity log 10 (^/cms —2 ) 

4.170(44) 

4.551(33) 

Absolute bolometric magnitude Mb 

4.09(7) mag 

6.09(7) mag 

Absolute visual magnitude My 
Separation between stars a/R© 
Distance d /pc 

4.19(8) mag 

31.06(3) 

4260(290) 

6.53(8) mag 


7 per cent larger while the secondary star is estimated 4 per 
cent smaller than in their study. These differences are most 
probably due to the higher accuracy of Kepler data relative 
to ground-based observations. 




0.48 0.49 0.50 0.51 0.52 0.98 0.99 1.00 1.01 1.02 

Phase Phase 



Figure 3. (a) Observed (dots) and computed (line) light curve 
of V568 Lyr. Primary (b) and secondary (c) minima expanded 
to emphasize the agreement, (d) Observed and computed radial 
velocities. 


4 MODELLING AND DISCUSSION 

We have simultaneously solved the radial velocity and light 
curve of the detached and relatively long period binary 
V568 Lyr and determined its orbital and physical parame¬ 
ters (Table 3). The well-detached configuration of V568 Lyr 
allows us to determine very accurate orbital and physical 
parameters (Tables 2 and 3). The distance to the clus¬ 
ter is found to be 4.26±0.29 kpc, equivalent to a distance 
modulus of 13.14 ± 0.22 mag. This is roughly consistent 
with other distance estimates. Harris & Canterna (1981) 
found a distance modulus of 14.0 ± 0.2 by analysing a 
CMD. Anthony-Twarog& Twarog (1985) found 12.6 using 
isochrones. Carraro etal. (2006) obtained 12.6 13.6 again 
with isochrones. Based on the fit to the white dwarf cool¬ 
ing sequence Bedin et al. (2008b) obtained distance scale 
as 13.50. Brogaard et al. (2011) determined 13.51 ± 0.06 
with binary system parameters. Basu et al. (2011) measured 
13.11 ± 0.06 for the red giant stars. 

To model this system and compare it with observa¬ 
tions allows us to test our understanding of stellar evolu¬ 
tion. V568 Lyr, in particular, gives us the opportunity to 
estimate the age of the cluster. We made with the ev code 
(Eggleton & Kiseleva-Eggleton 2002) and its much more 
powerful twin variant (Yakut & Eggleton 2005; Eggleton 
2006, 2010), both of which are based on the Cambridge 
stars code (Eggleton 1971, 1972, 1973; Pols et al. 1995). 

In single-star evolution, the effects of rotation and dy¬ 
namo activity on mass-loss and the proximity of the com¬ 
panion are often not considered. The ev code admits various 
non-conservative processes to be applied to the primary com¬ 
ponent of the binary system. In the newer twin variant, both 
components are solved simultaneously so that the effect of 
tidal dissipation on stellar rotation, orbital period and eccen¬ 
tricity and on dynamo activity, and hence on mass-loss, in 
both components are treated self-consistently. Stellar winds 
carry off angular momentum by way of magnetic braking 
and this angular momentum can be extracted from the or¬ 
bit if tidal dissipation is strong enough. We might imagine 
that tidal dissipation has played a role in the evolution of 


MNRAS 000, 1-6 (2015) 



























A Detached Binary Star in the NGC 6791 5 


V568 Lyr because its orbit is circular and many systems 
with P > 5d are markedly eccentric. Tidal dissipation leads 
to spin-orbit synchronization (or pseudo-synchronization in 
eccentric orbits) much more quickly than it causes circular¬ 
ization. 

It is not easy to determine the initial state of a binary 
system if the evolution is not conservative. We can reason¬ 
ably assume that the initial masses were larger than now 
but it is not clear, without some trial and error, by how 
much they were larger. Further, it is entirely likely that there 
is some non-linearity in the behaviour so that a straight¬ 
forward convergence procedure that assumes near-linearity 
may not in fact converge or at least converge badly. 

The metallicity of NGC 6791 has been estimated to be 
about two and a half times solar (Grundahl et al. 2008), 
rather surprisingly because the cluster is also estimated 
by them to be nearly twice the age of the Sun. How¬ 
ever the metallicity of the Sun has itself been subject 
to recent downward revision (Asplund, Grevesse, & Sauval 
2005) although there is probably still some uncertainty be¬ 
cause the new value seems to spoil the apparent agree¬ 
ment of earlier models of the Sun with helioseismology 
(Bahcall, Serenelli, & Basu 2005). We have used the tabular 
data of the OPAL data base (Rogers, Swenson, & Iglesias 
1996) for Z = 0.04, corresponding to nearly three times the 
value that Asplund, Grevesse, & Sauval (2005) obtained for 
the Sun, and to about 2.5 times the solar metallicity as es¬ 
timated by more recently by Caffau et al. (2011). However 
we adopt this value mainly because it leads to much better 
agreement between the observational and theoretical lumi¬ 
nosities, radii and masses of the components, as illustrated 
in Fig. 4. The zero-age helium abundance was taken to be 
Y = 0.32, and the mixing-length ratio to be 2.0. Aspects of 
the code have been described in detail by Eggleton (2006) 
and papers referred to therein. The only (slight) novelty is 
that convective core overshooting has been calibrated fairly 
carefully in a study by Eggleton & Griffin (2015) of 53 bi¬ 
naries almost all of which are evolved to giant or supergiant 
dimensions, where the issue of convective overshooting is 
much more challenged than by near-main-sequence stars in¬ 
cluding the present subject. The new model of overshooting 
in fact makes very little difference in the present case. 

After some experimentation, we evolved a pair of stars 
with initial masses 1.0994 and 0.8603 M0, spin periods of 2 
d each, an eccentricity of 0.3 and an orbital period of 15.79 
days. We expect the orbital period to decrease if the orbit 
circularizes, but it is also influenced by magnetic braking and 
by mass loss. In addition, it will also increase as the orbit ac¬ 
quires angular momentum from the spins of the stars by tidal 
friction. The particular model of the non-conservative pro¬ 
cesses in the code (Eggleton & Kiseleva-Eggleton (2002), 
but see below) led the system to evolve, after 7.733 Gyr to 
roughly the observed period, masses, radii and temperatures 
(cf Table 3): 14.42d, 1.0888 and 0.8285 M&, 1.428 and 0.781 
RQ, 1.949 and 0.309 L©, and temperatures 5712 and 4872 K. 
We see that the secondary component model is rather too 
hot, large and luminous, but probably not unacceptably so 
(2 - 3%). Fig. 4 shows the mean observed values (log L in so¬ 
lar units and effective log T in kelvin) as black squares, with 
an ‘error cloud’ of dots distributed according to an assumed 
Gaussian distribution with the standard errors as tabulated 
above. The blue circle is a point on the evolutionary track of 



3.75 3.7 

log T 


Figure 4. HRD of an evolutionary model with the TWIN code. 
Red, primary; green secondary. Black squares: mean observed val¬ 
ues. Dots: scatter due to random measurement errors; see text. 
Blue circle: a well-fitting model value of the primary. Blue aster¬ 
isk: co-eval point for secondary. 


the primary that is reasonably close to the observed value, 
and the blue asterisk is the coaeval point on the secondary’s 
(very short) track, i.e. a point on the same isochrone as the 
primary. However, in order to produce a nearly circular orbit 
we had to increase the model tidal friction by a factor of 100, 
and even then the eccentricity, starting from a hypothetical 
0.3, was only reduced by a factor of 6, which might just be 
compatible with the observations. We believe that the inad¬ 
equacy of tidal friction is likely to be a problem for any other 
model, because 14.5 d is a rather long period for two near- 
main-sequence stars to have circularized their orbit, even in 
quite an old cluster (Eggleton 2006). 

The fact that the two masses decreased by as much as 
1 per cent and 3 per cent in the course of 8 Gyr is perhaps 
noteworthy. Single stars similar to either component would 
lose less mass but that is because they would have spun 
down to much slower rotation rates by magnetic braking 
in the absence of tidal friction. The same non-conservative 
model as used above, when applied to the Sun, gives the 
observed solar-wind strength and Alfven radius, by design, 
at age 4.57 Gyr. This model predicts a loss of 1.4 per cent 
of the initial mass but almost all of this takes place in the 
first Gyr when the Sun is rotating in about 14 d or less. The 
lower-mass component of our binary loses mass faster, partly 
because its deeper convective layer is assumed to make it 
more active and partly because old K dwarfs are normally 
rotating even more slowly than the Sun. This one is rotating 
substantially faster if it is synchronized, which it surely must 
be. 
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